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ABSTRACT 

We discuss the low-curvature regime of time-dependent matrix theories proposed to 
describe non-perturbative quantum gravity in asymptotically plane-wave space-times. The 
emergence of near-classical space-time in this limit turns out to be closely linked to the 
adiabaticity of the matrix theory evolution. Supersymmetry restoration at low curvatures, 
which is crucial for the usual space-time interpretation of matrix theories, becomes an 
obvious feature of the adiabatic regime. 



1 Introduction 



Time-dependent matrix theories pQ have been introduced as an analogue of the Banks- 
Fischler-Shenker-Susskind flat space-time matrix theory [2] and of matrix string theory 
[31 HI [5] in an attempt to describe non-perturbative quantum gravity in time-dependent, 
possibly highly curved (or even singular) space-times. The original set-up of pQ has been 
later extended in various directions [7J EJ [91 [101 [TTJ [121 1131 El HI US]; in particular, a 
systematic generalization of the analysis to more general singular homogeneous plane-wave 
space-time backgrounds has appeared in [IB]. In close parallel to the flat space matrix 
theory conjectures, one may expect these models to give a complete quantum-gravitational 
theory of asymptotically plane-wave space-times. 

The usual construction of (time-independent) matrix theories [2] essentially relies on 
the type IIA superstring/M-theory duality conjecture. Namely, one compactifies a light- 
like dimension in the background 11-dimensional Minkowski space-time of M-theory, i.e., 
performs a discrete light-cone quantization (DLCQ) |17j . To make the construction more 
precise [18] (see also [IS]), one takes this compact dimension to be slighly space-like. A 
large boost relates this theory to M-theory with a manifestly space-like compactification on 
a very small circle. The latter system is related to weakly coupled type IIA string theory 
via the type IIA superstring/M-theory duality conjecture. Furthermore, for N units of 
momentum on the DLCQ circle and finite energies in the original reference frame used for 
DLCQ, N DO-branes have to be present in the type IIA theory, and the only surviving 
degrees of freedom are the massless open string states associated with the DO-branes. This 
set-up yields the matrix theory action. An analogue of this argument can be devised for 
the case of IIA superstring theory (rather than M-theory), yielding a matrix string action 
[21 SI E] (rather than matrix quantum mechanics) . 

Because plane waves enjoy a light-like isometry, the Minkowski space arguments can 
be generalized to plane-wave backgrounds, resulting in time- dependent matrix theories 
and matrix string theories. Even though, in the context of these models, novel physics is 
expected to emerge in the high-curvature regions of the time-dependent backgrounds, it is 
also important to understand in what precise manner the dynamics approaches classical 
space-time when the curvatures become small. 

A related issue is supersymmetry restoration. Supersymmetry is essential in the flat 
space matrix theory to protect the free propagation of gravitons, which in turn underlies the 
conventional space-time interpretation. In time-dependent matrix theories, supersymmetry 
is broken completely (whichever part is not broken by the plane wave backround, will 
be broken by the light-like momentum on the DLCQ circle). It is therefore crucial to 
understand why this does not prevent space-time from forming (at least in the low-curvature 
regime), an issue that has been raised since the original formulation of this class of models 
in [I]. 

Heuristically, the question of low-curvature dynamics has been addressed already in [TJ: 
it has been shown that, in an appropriate parametrization, the time-dependent matrix string 
theory reduces to a 2-dimensional supersymmetric gauge theory on a Milne space-time with 
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metric 

ds 2 = e 2r] {-drf + dx 2 ) , x ~ x + 2vr. (1) 

The supersymmetry is only broken by the identification x ~ x + 27r, and, since the radius 
of the Milne circle becomes large at late times (i.e., in the low curvature regime), one could 
expect that, for a wide range of processes, the supersymmetry breaking will become in- 
visible. In |20j, the effective potential for the matrix string variables has been computed 
in the weak coupling expansion of the gauge theory. Even though the computation is not 
technically valid at late times, formal extrapolation of the resulting expressions suggested 
that the effective potential decays at late times, a feature that could be indicative of su- 
persymmetry restoration. Another attempt to study late time (low background curvature) 
dynamics of the time-dependent matrix theories has been recently undertaken in |21j . 

Our present objective is to re-address these issues in a maximally clear and simple 
fashion. For the 11-dimensional case of [6], the matrix action previously presented in the 
literature shows steep time dependences at late times. This makes the emergence of a 
near-classical space-time quite puzzling, as it superficially suggests a strong explicit su- 
persymmetry breaking, among other things. We shall show that, treated in appropriate 
variables, the relevant time-dependent matrix theory approaches its flat space counterpart 
at late times (low curvatures). The manner of convergence is somewhat subtle, but, since a 
bound on deviations from the flat space theories will be given, the issues of space-time inter- 
pretation and supersymmetry restoration are automatically resolved. One can understand 
this situation in a different way: the superficially steep time dependences in the original 
matrix theory action actually enter an adiabatic0 regime at late times, which again connects 
the dynamics to that of a time-independent matrix theory. For the 11-dimensional case, 
this is simply an equivalent and less straightforward way to view the geometry-inspired 
variable redefinition that eliminates the time-dependences at late times. 

For the 10- dimensional case of [TJ [TB], however, there appears to be no canonical vari- 
able redefinition that eliminates the time dependence in the Lagrangian at late times. (For 
instance, the 10-dimensional theory of [TJ can be seen as the 11-dimensional matrix theory 
described in the previous passage, with an additional compactification. However, if we 
attempt to perform the same variable redefinition as in the 11-dimensional case, the com- 
pactification radius becomes time-dependent, making the theory awkward to study.) Still, 
one can establish the onset of an adiabatic regime in the late-time dynamics of these theo- 
ries, which again connects them to their flat-space counterparts. (In justifying adiabaticity 
for this case, we shall rely on the standard flat-space matrix string theory conjectures [5], 
which are essential to a meaningful interpretation of the time-dependent matrix theories in 
any case, and hence already implicitly assumed.) 

Note that in the present paper it is not our objective to trace the complete evolution 
of states from early to late times and show that a near-classical space-time emerges from a 
generic initial state (which we do not expect to be the case). Explaining the emergence of 
a near-classical space-time from initial conditions is a major problem in cosmology, which 

1 Adiabaticity has recently surfaced [22 in the context of quantum gravity in time-dependent back- 
grounds, though the precise setting differs substantially from ours. 
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we do not address here. Rather, we shall show that the dynamics of the relevant matrix 
theories approaches their flat-space counterparts, if studied at late times. This implies that, 
given a state with a late-time near-classical space-time interpretation, matrix theory can 
consistently describe its further evolution. Establishing this property is already non-trivial, 
and it is an essential pre-requisite for a more thorough treatment of cosmological scenarios 
in our framework. 

We shall start our exposition by analyzing the 11-dimensional time-dependent matrix 
theory, followed by the 10-dimensional matrix string theory. In these treatments, we shall 
perform the algebraic manipulations explicitly, leading up to the derivation of the necessary 
dynamical bounds. In the last section, we shall explain how our formal manipulations are 
related to quantum adiabatic theory. 

2 Matrix quantum mechanics 

We shall start by briefly reviewing the 11-dimensional (quantum- mechanical) matrix theo- 
ries introduced in [6] as simpler analogues of the matrix string theories of pQ. The relevant 
11-dimensional (M-theory) background has the form 

ds 2 = e 2ax+ (-2dx + dx- + (dx') 2 ) + e 2px+ (dx 11 ) 2 , (2) 

or, in terms of the light-like geodesic affine parameter r = e 2ax+ /2a, 

ds 2 = -2drdx~ + 2ar{dx i ) 2 + (2ar) /3/a (c/x 11 ) 2 . (3) 

This metric satisfies the 11-dimensional supergravity equations of motion if the constants a 
and (3 are related as = —2a, or = 4a. The fact that these relations need to be imposed 
will not be relevant for our present considerations (it is essential, however, for the general 
consistency of the corresponding matrix theories) . 

Since translations in x~ form an isometry of the above background, the usual DLCQ 
argument (proposed in [18] and adapted to the time- dependent case in [1]) can be applied. 
The result |6J is a matrix theory that can be expected to describe non-perturbative quantum 
gravity in space-times asymptotic to (13"1). The bosonic and fermionic parts of the matrix 
theory action, respectively, have the following form: 

Sb = J drTr {^(^) 2 + ^^OVf 11 ) 2 - j(2ar) 2 [X\ X>} 2 

-|(2ar) 1+ ^[X\X n ] 2 j, ^ 

S F = Jdr {i6 T D T 6 - Ry/2a^e T li {X i , 9] - R{2arf^ a 9 T ln [X n , 6}} . 

The problems we intend to discuss can be understood already at the level of the action 
01]). The space-time backgrounds implicit in (j4]) are supposed to feature a light-like singu- 
larity at r = and become progressively more classical at large r. Yet, the explicit time 
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dependences in (jl]) superficially become more steep, if anything, at large r. This issue cer- 
tainly deserves further clarification. Additionally, there are no supersymmetries explicit in 
fll]). Since supersymmetries are crucial for the free propagation of well-separated gravitons 
(and hence, a robust geometrical interpretation) in the flat space matrix theory, one should 
attempt to find an analogue of supersymmetry in (jl]) that would enforce a similar type of 
dynamics. 

To address these important issues, we first note that the metric ([3]) describes a plane 
wave and, with the coordinate transformation 

u = r, z l = V2^x\ z 11 = {2arf' 2a x 1 \ v = x~ + + ^ Z ^\ (5) 

Aar 

it can be brought to the Brinkmann form 

ds 2 = -2du dv — j- — .„ JK — -du 2 + (dz 1 ) 2 + (dz 11 ) 2 . (6) 

(2au) 2 

This parametrization forces the metric to manifestly approach Minkowski space-time for 
large values of the light-cone time, which strongly suggests that the large time dynamics 
of the corresponding matrix theory will likewise approach the flat space matrix theory, if 
treated in appropriate variables. (As we shall see below, the convergence towards this limit 
is somewhat subtle, but the naive expectation will prove well-grounded.) 

The matrix theory corresponding to ([6]) can be constructed from 'scratch', as the metric 
enjoys the same ^-translation isometry as ([3]). More straightforwardly, one can apply the 
matrix analogu^l of the transformations (jSJ) directly to (jl]) to obtain 

S B = j drTi [(D T Z 1 ) 2 + (D T Z U ) 2 ] - | [\Z\ Z>} 2 + 2\Z\ Z 11 ] 2 ] 

a 2 {Z l ) 2 -{P 2 -2am ll Y \ (7) 
(2ar) 2 J ' 

S F = jdr {i6 T D T 6 - R6 T ll [Z\ 6} - RQ T lxi \Z xl , 6}} . 

Note that the Brinkmann form of the matrix action has previously appeared [16] in the 
literature (for the 10-dimensional case). However, the action in [16] corresponds to bringing 
the 10-dimensional metric to the Brinkmann form, not the 11-dimensional metric (as we 
have done presently). Both metrics are of the plane- wave form. 

The action (J7]) only differs from the flat space matrix theory by a term decaying as 1/t 2 , 
thus one may expect that the late time dynamics will be approximated by the flat space 
matrix theory and admit the usual space-time interpretation. However, the decay is quite 
slow and one might be worried about whether it is sufficient to ensure convergence. 



2 Note that there is no matrix variable corresponding to v, since v becomes the DLCQ circle in the 
standard formulation of the matrix theory. The remaining transformations are linear, so the issue of 
matrix multiplication ordering never arises. 
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To illustrate these worries, one may look at the straightforward example of a harmonic 
oscillator whose frequency depends on time as 1/t 2 : 



x + —x = 0. 
t 2 



The two independent solutions to this equation can be given as t a and t 1 ~ a , where a is a 
fc-dependent number. These two solutions are obviously quite different from a free particle 
trajectory, even though the equation of motion approaches that of a free particle at late 
times. The reason for this discrepancy is the slow rate of decay of the second term in (jSJ). 

However, in a physical setting, one is only able to perform finite time experiments. That 
is, one has to specify the initial values x(to) = Xq, x(t ) = t> and examine the corresponding 
solution between t and t + T. The solution is given by 



x(t) 



x (l 



a 



v t 



1 - 2a 



t_\ Voto_ 
to) 1- 



XqCI 



2a 



l-a 



(9) 



One can then see that x(t + T) = xq + VqT + O(T/t ), i.e., it is approximated by a free 
motion arbitrarily well if the experiment starts sufficiently late. 

It may be legitimately expected that the finite time behavior of the full time-dependent 
matrix theory given by ((7|) will be approximated arbitrarily well by the flat space matrix 
theory at late times, just as in the above harmonic oscillator example. We shall now 
prove it by constructing an elementary bound on dynamical deviations due to a small 
time-dependent term in the Schrodinger equation. 

We start with the following Schrodinger equation: 



ij t \^) = (H + f(t)H 1 )\^) 



(10) 



where H and H% are time- independent, and rewrite it in the interaction picture (with 
respect to H ): 



|$) = e' iHo{t - to) |0 , 



We then proceed to consider 



4» 



f( t yH (t-t ) Hie -iH (t-t ) |^ 



d_ 

dt 



\m) - mo)) 



dt 



(12) 



= -*/(*) ((^(^0)1 e iHo( '-' o) i/ie-^ o( *-* o) |e(t)) - c.c.) . 
Integrating this expression between to and to + T and making use of standard inequalities 
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for absolute values and scalar products, we obtain: 

t +T 




e iH (t-to) Hie -iH (*-*o)^( to ) 



C.C.) 



t +T 



<2 / ^|/(i)| V /(| e ^o( t - t o)iy 1 e-^o(t-to) e (^))) 2 v /(| eW ))2 



(13) 



T 








Now, assume that f{t) approaches at large times and consider a fixed |£(to)) = |£q) ( so 
we consider the evolution with fixed duration T of the same initial state |£o) starting at 
different initial times to). In this case, the first factor in the last line becomes arbitrarily 
small for large t , whereas the second factor does not depend on t . We then conclude that, 
for sufficiently late times, the finite time evolution of the state vector will be approximated 
arbitrarily well by |£(t)) = const, i.e., by the evolution with f(t) set identically to 0. 

It is then a simple corollary of the above bound that the time-dependent matrix theory 
dynamics becomes approximated arbitrarily well at late times by the flat space matrix 
theory, and, in particular, the supersymmetry is asymptotically restored (with all the usual 
consequences, such as protection of the flat directions of the commutator potential, and 
free graviton propagation). 

A note may be in order: even though ( fl3l) does show that for any given experiment 
(fixed |£(£o)) = |£o) and fixed T) deviations from the flat space matrix theory become 
arbitrarily small for sufficiently large t , this by no means implies that the convergence is 
uniform with respect to |£o), which affects the value of the second factor in the last line of 
([131) . This is as it should be: for any t there will be some experiments that will be able to 
discriminate between the time-dependent and flat cases (with a given precision), yet, such 
experiments will have to become more and more specialized (and less and less possible) at 
late times. 

3 Matrix string theory 

The 10-dimensional "Matrix Big Bang" matrix string theories [H [16] essentially differ from 
the 11-dimensional case we have considered above by the compactness of one of the space- 
time dimensions. This feature precludes a straighforward variable redefinition that would 
reduce the Lagrangian to a sum of time-independent terms and terms decaying at large times 
(because the compactification radius would become time-dependent in the new variables). 
We shall therefore need to resort directly to adiabaticity-inspired arguments to establish 
the emergence of near-classical space-time far away from the Matrix Big Bang singularity. 

The aim of the construction of matrix string theories J3j HI [5] is to develop a non- 
perturbative description of quantum gravity in 10 dimensions (with the perturbative limit 
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of this description given by the usual perturbative type IIA string theory). For the original 
time-dependent Matrix Big Bang matrix string theory of p], the 10-dimensional geometry 
is asymptotic to the linear dilaton configuration: 

ds 2 st = -2dy + dy- + (dy') 2 , 

<P = -Qy + . [ ' 



To construct the matrix string theory for the background (TBI) , one first lifts the back- 
ground (|14p to 11 dimensions via the usual conjecture of type IIA/M-theory correspondence. 
The resulting 11-dimensional space-time is 

ds 2 = e 2Q y + ' z (-2dy + dy~ + (dy*) 2 ) + e~ 4 ^ / 3 {dyf , (15) 

where y is a coordinate along the M-theory circle. This is followed by the DLCQ compact- 
ification of the light-like f-coordinate, interpreted as the M-theory circle of an "auxiliary" 
type IIA string theory. A T-duality (23] then relates the resulting theory of DO-branes on 
a compact dimension, i.e., a BFSS-like matrix theory with a compactified dimension, to a 
more manageable theory of wrapped Dl-branes. This procedure has been carried out (in 
a slightly different but equivalent way) in [TJ and has been reviewed in [23] and [16]. The 
resulting matrix string action is 

s = /tr Qovn 2 +0 T ^+^^-^ , 

9ym 

with the Yang-Mills coupling gyn related to the worldsheet values of the dilaton: 

e -0(y+(r)) e Qr 

9ym = ^—j = — • (17) 

2nl s g s 2nl s g s 

A generalization of this set-up has been proposed |T6j|. One can start with a 10- 
dimensional power-law plane wave: 

ds 2 st = -2dy + dy- + g tl (y + Wdyi = ~2dy + dy- + $> + ) 2m *(c^) 2 

i 

= -2dzUz- + J2 ma( ( 7;~ V ) 2 (^ + ) 2 + E(^ a ) 2 ' (18) 

e 2< ^ = (y+) 3b /( b + 1 ) = ( z + ) 3b/{b+1) . 

Here, the first and the second line represent the Rosen and the Brinkmann form of the same 
plane wave, respectively. In order for the supergravity equations of motion to be satisfied, 
one needs to impose [16] 

^k-ih-J^. (19) 

i 

The original background of [1J can be seen as a b — > — 1 limit of the above space-time 
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The 11-dimensional space-time corresponding to (j!8p is 



ds 2 n = -2dudv + J2 u2ni ( d y 1 ) 2 + u 2 \dyf 

(20) 



2dudw + Y na{n \ l \ x a )\du) 2 + MLJl x 2 (rfM)2 + y (da . a)2 + (d 



with rij related to m 8 by 2mj = (2nj +&)/(&+ 1). Note that, since the Rosen and Brinkmann 
coordinates (first and second line in the above formula) are related by a ^-dependent 
rescaling of the transverse coordinates, the identification of the (compact) x-variable is 
M-dependent. Thus, even though the second line of (|20|) approaches a flat space-time at 
large values of u, the time-dependent identification of x makes an immediate application of 
the derivations of the previous section impossible. 

The usual formulation leads to a matrix string action, whose bosonic part is given, in 
the Rosen coordinates of f|T8|) . by 

S RC = j drdaTr f-lg-^^F^s - ^ gij (r)D a X%X^ 

+ \g 2 YM g lk (r)g jl (r)[X\X^][X k ,X l ] ) j . (21) 
with the transverse metric given by the first line of (118p and the Yang-Mills coupling by 

e -0(2/ + O)) r -3fc/2(6+l) 

g YM = — 7 - 1 = -j-, • (22) 

One can further transform this action to the Brinkmann coordinates of the original plane 
wave, given by the second line of f|T8l) . to obtain |16j : 



Sbc = j drdaTr (-^ 9 y 2 M F? a - ~ (D T Z a D T Z a - D a Z a D a Z a ) 



+\g 2 YM [Z\ Z h \\Z\ Z b ] + l -A ab {r)Z a Z b \ , (23) 

where = diag{m a (m a — l)}/r 2 . The latter form of the action only differs from a SYM 
gauge theory with a time-dependent coupling by the term involving A a b. The contribution 
of this term at late times can be shown to be negligible by an argument very similar to the 
one employed for the 11-dimensional case in the previous section. 

We thus end up at late times, both for the original Matrix Big Bang of pQ and for its 
generalization [16], with a super- Yang-Mills theory with a growing time-dependent coupling 
that is a power-law or exponential function of time, and we have to analyze this particular 
large time large coupling limit. Again, a puzzling feature here is that the time dependences 
remain steep at large values of r (set equal to y + by the gauge choice), far away from the 
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singularity of the original plane wave. The question of what happens at late times has been 
discussed in the literature [TJ |2TJ without definitive quantitative answers provided. 

The large coupling limit of the time-independent super- Yang-Mills theory is essential to 
the gravitational interpretation of the flat space matrix string theory [5J. Note that this 
limit is distinct from the one we have to consider (time-dependent versus time-independent 
theory), however, in the adiabatic regime that our considerations will establish, properties 
of the time- dependent theory are related to properties of time- independent "snapshots" of 
the time- dependent Hamiltonian. We shall now assume the standard conjectures about 
the large coupling behavior of the time-independent theory [5] and show that they result 
in an emergence of an adiabatic regime within the time-dependent theory. The relevant 
manipulations will be carried out in a pragmatic fashion, leading to a construction of the 
low-curvature limit. We shall explain the relation of these manipulations to the general 
adiabatic theory in the next section. 

In [5], time- independent super- Yang-Mills theories have been conjectured to converge 
in the infrared (gym — > oo) limit to a sigma model conformal field theory on a permutation 
orbifold (with target space S^R 8 = (R 8 ) N /S^, where Sn permutes the iV copies of R 8 ) at 
finite iV (and to second quantized free strings on Minkowski space in the large N limit). 
The considerations given in support of this claim are not rigorous, so we shall have to make 
natural assumptions below about which particular form of convergence is implied. 

Firstly, we assume that convergence occurs at the level of eigenstates: we postulate that 
there exists a family of eigenvectors \^f n (gyM)) of Hym{qym) suc h that 

+ £ ^ n ( gyM ^ = En (g YM ) \^ n (g YM )) ■ (25) 

k=1 9ym J 

Here \{A, 9(x)}) symbolically denotes eigenstates of the canonical coordinates of the super- 
Yang-Mills theory (the details depend on the quantization procedure). Their appearance 
in the above expressions is necessary to compare wave-vectors of super- Yang-Mills theo- 
ries at different values of the coupling (which are technically defined in different Hilbert 
spaces). Further, \]>p ermorb - an d _gp crm - orb - denote wave-functionals of energy eigenstates of the 
permutation orbifold conformal field theory and their corresponding energies. (The above 
set-up assumes that the spectrum at finite gyu is a continuous deformation of the limiting 
spectrum at gyu ~ > oo. This may not always be so, for example, a discrete spectrum of 
bound states at finite gyM may merge in the continuum at gyu — > oo. However, we are not 
aware of any such states in the context of perturbative string theory, nor do we anticipate 
that they could have any dramatic effect on our considerations. We shall therefore assume 
( 121)1) for our present purposes.) Note also that our summation over n is a symbolic nota- 
tion that implies integration over the continuous spectrum and summation over all discrete 
eigenvector labels, and none of our derivations use any assumptions about discreteness of 
the spectrum. 

If the Hilbert space were finite-dimensional, the above conditions would guarantee that 
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the dynamical behavior of the theory reaches the free string limit suggested in [5] at large 
values of Qym- However, for an infinite-dimensional space of states, additional conditions 
enforcing not-too-poor convergence of highly excited eigenstates need to be imposed. In- 
deed, in a physical setting, in order to declare that one theory approaches another in a 
certain limit, one needs to ascertain that the evolutions of the two theories become arbi- 
trarily close to each other for any finite energy normalizable initial state and for any finite 
duration of the experiment. Furthermore, this convergence needs to be uniform for all nor- 
malizable initial states with energies below some fixed value and all experiments lasting less 
than some given duration (for, were it not so, there would always exist some finite energy 
finite duration experiments discriminating between the two theories with a finite precision). 
We embody these conditions in a statement that 

oo , 

J2 c n e- iE ^ T ({A, 9{x)}\M9ym)) = £ c n e^ E ~ h ■ + £ (26) 

n n k=l Sym 

with ^2 \c n \ 2 = 1 and ipk being a set of normalizable wave-functionals. We further impose 
that the norms of ipk are uniformly bounded for any T less than a chosen fixed value and 
for any initial state energy ^ E , £ erm ' orb '|c n | 2 less than a chosen fixed value. 

We shall assume something slightly stronger than ( |26i) . Namely, we shall assume that 
the Yang-Mills coupling constant in the energies and eigenvectors on the left-hand side of 
(1261) can be sent to infinity independently, still resulting in a power series expansion: 

oo . 

£c n e-^) r ({AW}^ £ Tr, (27) 

n n k,l=o -J^^z 

with the same type of uniformity specifications we made under f )26|) . We shall further 
assume that (|27|) can be differentiated with respect to T without losing uniformity, which 
yields (with T set to and g 2 sent to infinity): 

oo , 

c n {E n {g YM ) - E P~ h -)W~ h - = ( 28 ) 

fc=i 9ym 



As before, all these conditions are identical to (125]) for any finite-dimensional space of states, 
but in an infinite-dimensional space, they constrain the convergence (in Qym) a t large n. 
(These conditions are still considerably milder than demanding uniformity of convergence 
in n, which would be quite unphysical.) There is a possibility that the derivations below 
can be made without relying on anything beyond ( 1261) . but we have not been able to devise 
such an argument. 

Now we turn to the Schrodinger equation of the time-dependent super- Yang-Mills theory 

ij t |*) = H YM (g YM (t)) |tf> (29) 
and expand the state vector in the basis given by (|25|) : 

l*) = £^(*)l*»(s™(*))>- (30) 
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This yields 



i^jr + % ^ c m (t)g Y M (^n(gYM)l j^— \^ m (9YM)) = E n (g YM )c n (t). (31) 

We can rewrite this equation as 

dc 

i~£ = E P n — h -C n + HnmCm(t), (32) 

m 

where H nm are bounded by a negative power of gyu'- 

H nm = 0{l/g YM ), (33) 

with the precise value of the power depending on whether gyu is a power-law or exponential 
function of t. Note that the relevant coupling dependence coming from the second term on 
the left-hand side of fl3T|) is gYAi/gyM^ whereas the energy corrections on the right-hand 
side of fl3Tj) are of order 1/gyM- 

The fact that the matrix elements H mn are decreasing functions of time suggests (though 
by no means in a conclusive way) that the evolution approaches that of the permutation 
orbifold CFT at late times. We shall now analyze (13"2"1) in more detail to establish that this 
convergence does indeed take place. 

Consider first a Schrodinger equation of the form 

<^|*> = (flo + J7i(t))|*>, (34) 

where H is time-independentH. In parallel to the derivations of the previous section, we 
can rewrite it in the interaction picture (with respect to H ): 

|$) = e - iHo{t - to) \0 , i^- \0 = e iHo(t~to) Hi ^ e -iH (t-t ) |^ _ ( 35 ) 

We then proceed to consider 



d 
dt 



!£(*)> - KftO) 2 = -^((^o)l^)) + c.c.) 

ai (36) 



= -i (^{t^e^-^H^e-™^-^ \£{t)) - c.c.) . 
Integrating this expression between to and to + T and making use of standard inequalities 



3 It is an interesting mathematical question, to which the authors do not know the answer, what precise 
conditions should be imposed on Hi(t) in order to make its contribution to the evolution small. There 
are many notions of convergence for Hilbert space operators, and saying that Hi(t) becomes small in some 
limit is vacuous, unless the precise manner of convergence is specified. 
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for absolute values and scalar products, we obtain: 



t +T 

|£(i + T)> - |e(*o)> f = ~i J dt ((e^-^H^e-^-^ato) |£(f)> - c.c.) 

to+T 

< 2 / tfti 



^lle^o^O^^e^o^o) l^to))!! 2 ^!! |£(t)>|| 2 (37) 



to 



2 J dtsJ\\H x {t + t )e- iH ^\^{t Q ) 



For the Schrodinger equation (l32|) . the state vector is given by the set of numbers {c n }, 
and we take H to be the first term on the right-hand-side, and H\(t) to be the second 
term (the detailed expression for H\ can be read off from (I3lj) ). We furthermore denote 
|£(£o)) = { c n°' ) }, and observe that 



i/i(t + to)e- iffo *|e(to)) 



= J" (qym (* n (ffYM{t + *o))| -7— — \y m (9YM(t + t ))) exp [-lE^^t] c$ ) 

™ \ d 9YM J 



(E n (g YM (t + t )) - ) exp [-z^«»-« b -t] . 



The norm of the first term (i.e., the sum over n of the square of its absolute value) is equal 
to the norm of the Hilbert space vector 

9ym V] eg? exp [-iE^H] -1— \V m (9YM(t + to))) (39) 
m d 9 YM 

(we use ^2n\^n(9YM{t + to))} (^n{9YM{t + 1 ))\ = 1) and it is bounded (uniformly with 
respect to t) by the g 2 -derivative of the condition (j2"T|) with gi — > 00 to be 0(gyM/ 'Oym)- 
The norm of the second term is bounded (uniformly with respect to t) by ( 128]) to be 
0(1/ gYhi)- Since the bounds are uniform with respect to t, they can be immediately 
integrated in ( 157)) . 

We then conclude that 

\at + T)> - |£(t )> f = TO(l/fi M (to)), (40) 

i.e., deviations from the free string on the permutation orbifold become arbitrarily small at 
large to- Just as in the previous section, supersymmetry restoration at late times becomes 
a simple corollary of the above bound. 
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4 Adiabaticity 



In the preceding sections we have derived bounds controlling the convergence towards 
the late-time limit of time-dependent matrix theories. A reader familiar with quantum 
adiabatic theory (see, e.g., [25]) will immediately recognize the structure of our manip- 
ulations. Indeed, for both 11-dimensional and 10-dimensional cases, we have related the 
time-dependent theory to properties of time-independent theories (flat space matrix theory, 
super- Yang-Mills theories), which is characteristic of the adiabatic approximation. In this 
section, we shall briefly review quantum adiabatic theory and display its connections to the 
derivations in the previous sections. 

Given a general quantum system with a time-dependent Hamiltonian H(t) and a state 
vector satisfying the Schrodinger equation 

i^|*> = J ff(t)|*>, (41) 

one can always expand the state vector in a (time-dependent) basis of instantaneous eigen- 
vectors \^n(t)) of H(t): 

l*> = $> n (t) l*n(t)> , H(t) |* n (t)) = E n (t) |tt n (t)) . (42) 

n 

Note that our summation over n is a symbolic notation that implies integration over the 
continuous spectrum and summation over all discrete eigenvector labels, and none of our 
derivations use any assumptions about discreteness of the spectrum. The Schrodinger 
equation then takes the form 

^ + J t l*™(*)> = E nit)Cn{t). (43) 

m 

An adiabatic regime occurs when the second term on the left hand side becomes small. 
Heuristically, this happens when Hit) varies slowly in some sense, and so do |\I/ n (t)), so 
that their derivatives can be neglected. In general, it is difficult to spell out more handy 
conditions for the emergence of this regime. However, in particular cases, simple and explicit 
adiabatic parameters can be constructed, as we shall see below. 

If the second term on the left hand side of (j4"3l can indeed be neglected, the equations 
can be solved trivially to yield 



c n (t) = C n exp 



J dtE n {t) , |^) = ^C„exp -i J 



dtE n {t) 



!*»(*)> • (44) 



Note the close similarity between these approximate adiabatic solutions to the time-dependent 
Schrodinger equation and the familiar solutions to a time-independent Schrodinger equa- 
tion: \^) = 5^ n C n exp [— iE n t] l^n)- (The stationary eigenvectors |\& n ) are simply replaced 
by the instantaneous eigenvectors \^ n (t)), and the E n t in the phase factors are simply 
replaced by J dtE n {t).) 
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The relation between the time-dependent and time-independent systems becomes even 
more straightforward if the instantaneous spectrum E n (t) of H(t) scales uniformly as a 
function of time, namely, 

E n (t) = \(t)E$>\ (45) 

where the E^ do not depend on time. In that case, the phase factors in ( 1441) become simply 

—iEn^ J dt\(t) , in other words, they differ from the stationary case only by replacing t 

with J dt\(t), independently of which state vector one is dealing with. 

For the class of systems characterized by (|45|) . it is convenient to perform a variable 
redefinition in the Schrodinger equation. (145]) implies that there exists a time- dependent 
unitary transformation S(t) such that 

&(t)H(t)S(t) = X(t)H , (46) 

where Hq does not depend on time and possesses the spectrum E^ ■ One can then introduce 
|$) = S^(t) |^), satisfying 



ij t m+i(s\t)p(t)\m = x(t)H \^. 



(47) 



An adiabatic regime occurs when the second term on the left hand side can be neglected 
compared to the right hand side. The time dependence responsible for this relation between 
the two terms can be isolated into the operator 

S\t)±S(t)) . (48) 



X(t) V dt 

If the second term on the left hand side of f|4T|) can indeed be neglected, the equation is 
solved to yield 



= S(t) |$(t)) = S(t)exp 



t 



-iH J dtX{t) 



to 



S\t ) |*(t )> (49) 



The evolution is essentially that of a stationary system described by Hq, except that a 
unitary transformation S(t) is performed and the "time flow" is deformed from t to J dtX(t). 

The situation simplifies further if S(t) corresponds to a (time-dependent) linear trans- 
formation of the canonical coordinates 

Qk = s k i(t)qi (50) 
(and the corresponding transformation of the canonical momenta): 



S(t) = Vdet s / dq \sq) (q\ = = f dq \q) (s 1 q\ . 
J v det s J 



(51) 
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Differentiating S(t) with respect to t is somewhat subtle and is most conveniently performed 
using 



Then, 



= i(s x ss 1 q) k p k S(x - s x g) = i(s 1 s) kl p k q l |s x g) . 

^5 = 5 - iCs-^W*) = 5 Li( s -is) w { ft , Pfe }^ . (53) 

Hence, (|4"T|) takes the form 

4 1$) + -(s^s)^,,?*} |$) = A(t)i7 1$) • (54) 
dt 2 

All the time dependences have now been isolated into numerical pre-factors of the operators. 
The characteristic ratios (adiabatic parameters) quantifying the neglect of the second term 
on the left hand side with respect to the right hand side can now be given9 as a matrix: 



X(t) 



(55) 



One can examine how the above formulas work for the familiar case of a time-dependent 
harmonic oscillator with H(t) = {p 2 + u 2 (t)x 2 )/2. If one makes use of (1511) corresponding 
to the transformation x = x/y/u, one obtains Sft(t)H(t)S(t) = u{p 2 + x 2 )/2. In other 
words, our above analysis applies with s = 1/a/cJ and A = u, so that fl55|) becomes simply 
proportional to u/u 2 , i.e., the relative change of u per period of the oscillations, which is 
the familiar adiabatic parameter of the harmonic oscillator. The new Hamiltonian after the 
transformation has been effected can be read off (|54[) as 



^ , 2 , 2-s ^ 



H = —(p + x ) + — (px + xp). (56) 

We now turn to the derivations of the preceding sections. The variable redefinition 
([5]) we have employed for the 11-dimensional matrix theory is precisely of the form ( )50|) . 
and it converts the time-dependent system to a time-independent system plus a correction 
decaying at large times. If we forget about the geometrical interpretation of this variable 
redefinition, it simply becomes a particular case of ( 1501) demonstrating that the (superficially 
steep) time dependences in the Rosen form matrix theory action (j3J) in fact become adiabatic 
at late times, and the system is well approximated by its time-independent counterpart, 
i.e., the flat space matrix theory. 

For the 10-dimensional case, equations ( 13011311 are precisely of the form ( 1421143]) . Fur- 
thermore, ( T45|) is almost satisfied at late times (because the spectrum approaches a constant 



4 For the simple systems we are now considering, a purely classical consideration involving the trans- 
formation q k — would produce the same adiabatic parameters. The way we have derived them 
here implies automatically their validity for the quantum case, which is the physically relevant regime for 
gravitational matrix models. 
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limit). Our analysis of the 10- dimensional case shows explicitly that the second term on the 
left-hand side of fl3Tl) can be neglected in comparison to the other terms (provided that the 
convergence of time-independent matrix string theories to the Dijkgraaf-Verlinde-Verlinde 
limit is sufficiently tame), which is by definition the adiabatic regime. 

One may try to object that adiabaticity is not a relevant term for our discussions, 
since a variable redefinition brings the equations of motion to the form where all the time 
dependent terms become small at large times. This objection is vacuous, however, since it 
would also apply to a wide range of systems commonly thought of as adiabatic. Indeed, 
if E n (t) in (143 p approaches a constant limit at late times, ( I43p will take the form where 
all the time dependent terms become small at large times whenever an adiabatic regime 
occurs at late times, in this general setting. Likewise, the familiar time-dependent harmonic 
oscillator H{t) = {p 2 + uj 2 {t)x 2 ) /2, the simplest system used for text book demonstrations 
on adiabaticity, can be converted to 

H{t) = V -±^- + (p X + xp) , (57) 

if we start from (|56|) and introduce a new time variable r = J dtu(t) (the dot in the 
above formula still denotes the derivative with respect to the old time t to maintain the 
familiar expression for the adiabatic parameter, u/u 2 = d T u/u). Then there is a one-to- 
one correspondence between the adiabatic regime and the (new) Hamiltonian being almost 
constant: both occur when the adiabatic parameter u/u 2 is small. Yet, this mathematical 
structure does not prevent anyone from employing the term 'adiabatic' for the adiabatic 
regime of a time-dependent harmonic oscillator. 

There appears to be widespread intuition that adiabaticity in quantum mechanics is 
somehow connected to discreteness of the energy spectrum, and adiabatic parameters 
emerge from comparisons of the rate of change of various terms in the Hamiltonian to 
energy spacings in the discrete spectrum. This intuition stems from the simplest versions 
of adiabatic theorems (proved, for example, in |25j), which are sufficient, but certainly 
not necessary conditions for adiabaticitjo To dispel the doubts regarding adiabaticity in 
systems with a continuous spectrum, one may simply notice that an inverted harmonic 
oscillator H(t) = (p 2 — u 2 (t)x 2 )/2 can be brought to the form 

2 2 

H(t) = V -Z^_ + £_( px + xp ) (58) 

by the same transformations we used in obtaining ( |57|) . This system will be adiabatic 
whenever Cj/uj 2 is small by virtue of the bound (1131) . irrespectively of the fact that the 

5 One must also keep in mind that there are different ways to reach the limit of slow relative variation of 
the Hamiltonian, and they may result in different adiabatic parameters. For example, in adiabatic theorems 
proved in [25] , one considers a general Hamiltonian H(sT) where s changes between and 1, and the limit 
of slow variation is reached by sending T to infinity. The evolution is examined between sT = and 
sT — T, and demanding adiabaticity is a very strong requirement, since the non-adiabatic terms should 
produce a negligible contribution even over the huge time interval T. We study adiabatic evolution on hnite 
time intervals for Hamiltonians of a very particular form. 
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spectrum is entirely continuous. For the matrix theories we have explored in this article, 
adiabaticity has been proved by constructing explicit bounds on deviations from the strictly 
adiabatic evolution. 

5 Conclusions 

We have considered the low-curvature regime of time-dependent matrix theories and matrix 
string theories, and displayed the relation between the emergence of near-classical space- 
time and adiabaticity of the time dependences in the matrix theory actions. In this context, 
supersymmetry of the matrix theories (explicitly broken by the time dependence) is natu- 
rally restored at low curvatures, and the conventional space-time interpretation of matrix 
theories becomes viable. 
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